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Abstract 

We investigate theoretically polyelectrolyte bridging interactions on the two-body 
level. The model system is composed of two macroions with two oppositely charged 
flexible chains. The electrostatic interactions are treated on the Debye - Hiickel level. 
The formal level of the theory is provided by the Feynman-Kleinert variational method 
generalized to include also self-interactions between polyelectrolyte segments. The 
variational equations are shown to exhibit two solution branches corresponding to 
strong and weak coupling, whereas conformations of the chain can be described as 
weakly or strongly paired. We investigate the effective pair interaction between the 
macroions in the parameter space and comment on the relevance of the calculation 
for bridging interactions in experimental context. 

1 Introduction 

Polyelectrolytes are ubiquitous in colloidal systems and play a fundamental role in de- 
termining the interactions between as well as stability and structure of various molecular 
assemblies pQ. Their effect on colloidal interactions has been studied and exploited in 
various technological contexts ranging from paper industry to pharmaceutical industry 
0. It seems however that their most basic role is played in the biological context where 
their importance can be hardly overestimated. They are an essential and fundamental 
component of the cellular environment and make their mark in its every structural and 
functional aspect (Hj. 

The behavior of polyelectrolytes in biological context has without any doubt been one 
of the foci of the soft matter research for quite a few years now 0j . The intense work on the 
interactions and mesophase behavior of the most studied polyelectrolyte in the biological 
context if not in general, i.e. DNA, has elucidated many fascinating physical aspects 
of this molecule and the repercussions that they have on the structure and function of 
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biological matter [jj]. The mesoscopic interactions between many DNA molecules (Bj and 
elastic properties of single DNA molecules in different solution conditions [7j have been 
measured directly and are understood on a fundamental physical level. Not all biological 
polyelectrolytes or all solution conditions have been or indeed can be studied at quite 
the same level of detail. Sometimes a lot less information then direct measurement of 
molecular interactions at all macromolecular densities is experimentally available. Studies 
at low macromolecular densities in systems where polyelectrolyte behavior is expected to 
show its mark usually only lead to second virial coefficients and not complete interaction 
curves [HI El El- This is not due to poorly designed experimental setup but shows a rather 
fundamental limit in the amount of information that can be provided by experiments 
at these conditions. In this situation one has to rely heavily on different models of the 
mesoscopic interaction potential and the way they transpire through the measured second 
virial coefficient. This is the first point of departure of this work. 

The second one is a very peculiar interaction in polyelectrolye systems, where long 
charged polymers can mediate interactions between macroions of opposite charge 
(and references therein). The term bridging interactions is usually applied to this sit- 
uation where a single chain can adsorb to different macroions and via its connectivity 
and elasticity mediate attractive interactions between them. These interactions have been 
studied intensively both experimentally as well as theoretically I14| IT ^ I16>] . Surface 
force apparatus and atomic force microscopy have provided direct data on the separation 
dependence of the bridging interaction between macroscopic surfaces with polyelectrolyte 
chains either grafted or in chemical equilibrium with a bulk solution ^2] • Theoretical work 
has added a clear mesoscopic picture for the bridging interaction between macroscopic sur- 
faces and elucidated the effects of salt and non-electrostaic excluded volume effects on the 
strength and range of this interaction |13| I14| 115] 116] . Since it is based on sometimes severe 
model or formal restrictions there is no single theoretical approach that is able to account 
for all experimentally observed details or is able to explore in comparable details all the 
regions of the parameter phase space [11] . The fact that the effects of bridging interaction 
between small macroions ^3] as opposed to macroscopic surfaces have been studied 
to a much smaller extent thus makes a strong point for its reevaluation. 

The main motivation for this task are recent experiments on the second virial coefficient 
of the nucleosomal core particles (NCP) in low density regime at various solution conditions 
[101 117] . NCPs represent the lowest level of the chromatin organization in eucaryotes and 
have recently been resolved at an atomic resolution JH] • They consist of a histone protein 
octamer core with 146 bp of DNA tightly wrapped around it, giving it an approximate 
cylindrical shape of a radius ~ 55 A, a height of ~ 57 A and a structural charge of ~ -250. 
This complex is stable in an aqueous solution from mM to 750 mM monovalent salt ionic 
strength [17] . The charged histone N-termini or N-tails can desorb from this complex and 
basiaclly play the role of grafted flexible polyelectrolyte chains of an approximate total 
charge of ~ +90. They remain essentially adsorbed to the DNA segment of the NCP at low 
ionic strength but tend to assume a more extended conformation as the ionic strength is 
increased [10] . The application of classical an d manometric :2Q, osmometry provided 
data on the osmotic pressure of the NCP in NaCl solutions of variable ionic strength from 
which the second virial coefficient was deduced quite accurately ^U] . It was demonstrated 
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that the second virial coefficient is a non-monotonic function of the ionic strength and a 
plausible though tentative hypothesis was given, that the non-monotonicity of the second 
virial coefficient might be due to the bridging interaction mediated by the extended N-tails 
of the NCPs jlUj . Very similar non-monotonic dependence of the second virial coefficient 
was seen also in apoferritin solutions where a bridging interpretation is hard to envision 

Motivated by this important experimental result and its tentative interpretation we 
embarked on a detailed study of the interaction between charged macroions with grafted 
oppositely charged polyelectrolyte tails, as a function of the ionic strength of a monovalent 
bathing salt solution. The level of the theoretical calculation had to be considered, that 
would allow for a straightforward evaluation of the second virial coefficient, being an 
integral of the underlying interaction potential, and its dependence on the ionic strength. 
Also since the N-tails in the motivating experiment can not be considered as infinitely 
long, i.e. there are finite size effects that need to be taken properly into account, it seems 
that a mean-field theory of bridging interactions as formulated for the case of interacting 
charged planes and based on the ground-state dominance ansatz |13|ll6j can not be simply 
implemented to the present case. Finite size effects are quite difficult to deal with on the 
MF level, especially if one needs to evaluate the interaction potential between macroions 
in a very large (ideally infinite) range of separations. In view of all this we formulated 
a variational |2J , two particle (two macroions) theory of the bridging interaction that 
starts from an explicit mesoscopic ployelectrolyte model and includes the interactions of 
the polyelectrolyte chain with the macroions, the interactions of the chain with itself and 
connected with it the effect of the electrostatic stiffening of the chain, as well as the 
configurational entropy of the chain. In some respects the theory proposed here could be 
viewed upon as a variational version of the Asakura - Oosawa theory [22]. The finite size 
of the chain can be relatively straightforwardly dealt with on this level of the theory and 
gives rise to important features of the two-particle bridging interaction that are lost in the 
simplest, ground state dominance formulation of the MF theories. 

The organization of the paper is a s follows: we will first describe the model and give an 
introduction to a modified variational Feynman-Kleinert approach to the polyelectrolyte 
chains. We will derive the main equations and solve them numerically for different condi- 
tions. We will show that in general the bridging interactions for this model system comes 
in two varieties that we dub the strong and the weak coupling limit. The form of the 
total interaction between the macroions will be obtained numerically for all values of the 
intermacroion separation as a function of system parameters such as the amount of fixed 
charge on the macroions, the length of the polyelectrolyte chain and the screening length 
of the intervening bathing solution. We will discuss the ramifications of these results for 
the salt dependence of the second virial coefficient and point to the possible shortcom- 
ings of the calculation and guess about the way to possibly circumvent them. We will 
finally comment on the significance of the present calculation for the understanding of the 
bridging interaction in the NCP solution system. 
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2 Model 



The model system that we take as the starting point of our theoretical investigation of the 
bridging interaction is quite simple: it is composed of two spherical point macroions with 
M negative fixed charges plus two oppositely charged chains, each with N monomers, one 
per each monomer. The pair interaction potential u(r',r) between all the charges in the 
system will be taken of the screened Coulomb (Debye-Hiickel) form [Q 

g 2 e — re|r-r'| eg / 

u ( r '' r ) = i 7T or else u ( k ) = — 7H~, — 2T' W 

47reeo |r — r| eeo(k z + k 2 ) 

in real in in the Fourier space, the form we will need later on. k is the inverse Debye length, 
eo is the elementary charge, one per each Kuhn's length and the rest of the notation is 
standard. Obviously counterions are not explicitly included in this model. The interaction 
potential between the polyelectrolyte and the macroion charges is assumed of a similar 
form viz. 

„ „ _— k\t— ri| „ „ D — k\t— ral 

eie e 1 11 e 2 e e 1 ^ 

(t>ext{r) = -. i r + -. i r + • • • , (2) 

47reeo |r — ri| 47reeo |r — r2| 

where ri, r 2 etc. are the positions of macroions and their charges are e\ = e 2 = Meo etc.. 
Our model is thus a very straightforward generalization to many macroions of the model 
used in polyelectrolyte adsorption studies [23] . 

We will use a standard (Edwards) model |2l] for the polyelectrolyte chain where the 
mesoscopic Hamiltonian has contributions from chain connectivity, interactions between 
the segments of the chain and the interaction with an external field due to the presence 
of two macroions. It is written as 

3 2 r N -y 2 r N r N 2 r N 

/3H[ri(n)] = — g V) / rf(n)dn+-/3 V / / u(r -(n), rj(n))dndn'+(3 V / 4>ext(r%in))dn, 

(3) 

where ^ is the Kuhn's length, n(r', r) is the pair interaction potential , while 4> e xt( r ) is the 
external interaction potential. The indexes i,j stand for the two polyelectrolyte chains. 
Fig. Q schematically presents the mesoscopic model on which the present evaluation of the 
bridging interactions is based. Clearly the non-pairwise additive effects such as bridging 
between multiple macroions mediated by a single chain, have been completely disregarded 
in this model. The finite macroion size effects have also been disregarded. Also the model 
is based on a linear theory (Debye-Hiickel) of Coulomb interactions and thus can not 
capture non-linearities such as charge renormalization or counterion condensation. It can 
however take into account the electrostatic stiffening of the chain as well as the finite chain 
size effects. 

The grafting of the chains to their respective macroions is not taken into account 
explicitely on the Hamiltonian level. First of all in this model system the effects of grafting 
are small or indeed negligible |15j to the extent that they are always overpowered by the 
much stronger electrostatic interactions. This would of course not be the case for e.g. 
electrostatic brushes where grafting has to enter the description of the chain already at 
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the Hamiltonian level. The grafting of the chains is only taken into account via their 
center of mass coordinates in the way explained later. 

3 Variational ansatz and formalism 

Since the interactions with external fields as well as the self-interactions along the polymer 
chains are highly non-linear in terms of their spatial dependence, the statistical integral 
corresponding to Eq. is in general impossible to evaluate analytically and is difficult to 
evaluate even approximately. Instead of taking recourse to a numerical approximation we 
will rather introduce a harmonic variational ansatz |25| I26j that will make the evaluation 
of the statistical integral straightforward. The parameters of the variational ansatz will 
be chosen so as to minimize the upper bound of the exact free energy. This procedure 
is generally known as the Feynman - Kleinert variational method 26 : and its application 
to the case of self-interacting polymer chains in the bulk has been discussed in detail 
elsewhere [25- It was shown that it captures all the salient features of polyelectrolyte 
behavior in the bulk. 

The details of the variational approach are as follows [25 : the two polyelectrolyte 
chains are trated as Gaussian blobs positioned at their centers of mass, i"oi and ro2 respec- 
tively. The width of the Gaussian blobs is determined variationally from the strength of 
the external fields as well as self- interactions along the chain . The Gaussian blobs thus 
behave as effective "particles" with finite extensions. In this respect we can claim that 
the present theoretical framework represents a kind of Asakura-Oosawa theory where the 
effective size of the macroions is determined variationally. The final statistical integral is 
then obtained by integrating over the two centers of mass roi and ro2, i.e. the coordinates 
of the effective (Asakura-Oosawa) "particles", or indeed by finding the configuration of 
the centers of mass that gives the largest contribution to this integral |27| . 

For the variational ansatz corrsponding to two effective Gaussian chains we will chose 
a general harmonic hamiltonian of the form 

3 2 f N 3 2 f N 

0Ho[ri(n)] = r^V r?(n)dn + - V CjW) / (r<(n) - r 0i ) 2 dn + /3iVT(r i, r 02 ), 

21 i=i J ° 2 i=l J o 

(4) 

with periodic boundary conditions for Vi(n). This ansatz is obviously still dependent on 
roi for i = 1,2, that stand for the centers of mass of the two chains i.e. = j? f Ti(n)dn, 
as well as the functions Ci( r 0i) and £(roi,ro2) that will be determined variationally. 

The term with Cfi r 0i) obviously represents an external harmonic potential, centered 
on roi, that acts either to confine or expand the chain, depending on its sign. A simple 
limiting form of ( 2 can be derived only for the case of a single self-interacting chain and 
is given in [2j|. In Eq. 0]this term was taken with a positive sign but we will argue later 
that it can as well be negative. The term /?iVX(roi, 1*02) simply represents the value of this 
harmonic external potential at the centers of mass of both chains. Again a simple limiting 
form of (3NC can be derived only for the case of a single self-interacting chain [25- As 
will become clear when we proceed, both quantities depend in a complicated way on the 
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interactions between the monomers as well as on the interactions between the monomers 
and external macroions. 

The statistical integral for the variational ansatz can be obtained in the following form 

E (N) = J V[n(n)} J V[r 2 (n)] e -/»«oM»)] = J J d 3 r 01 d 3 r 02 e -^o(r 01 ,r 02 ^ (5) 

The two polymer chains are thus represented as two effective Gaussian "particles" with 
an effective Hamiltonian given by .Fo( r oi> r 02)- The details of the implementation of the 
Feynman - Kleinert ansatz [22j for (self) interacting polymer chains have been given before 
|21| and we will rely on the formal developments described in that work. First of all let 
us introduce the radius of gyration defined as 
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where C(x) = cothx — ^ is the standard Langevin function. One can then derive [2^ 
that minimization of the upper bound of the variational free energy with respect to the 
function £(roi,i"o2) leads to the following equation 

/3iV£(roi,r 02 ) = -f £ 2 = l (?(r 0i )Na*(r 0l ) + 0W(r Q1 , r 02 ), (7) 

with 

PW(r 01 ,r 02 )) = P J d 3 r 4> ext (r)p(r,r 01 ,r 02 )++^(3 J J d 3 rd 3 r'p(r,r 01 ,r 02 )u(r,r')p(r',r 01 ,r 02 ) 

(8) 

W(rQi,ro 2 ) obviously represents the total interaction free energy, due to self as well as 
interactions with external fields, of a smeared monomer cloud with a Gaussian density 
distribution. p(r, i"oi,ro2) stands for the combined monomer density function of the two 
chains and has the form 

p(r,r i,r 02 ) = p a 2 (r, r i) + p a 2 (r, r 02 ), (9) 
where for each of the chains the monomer density distribution function is given by 
_/\r |j- _ ro; I 2 fc 2 

Pa ^ r ' roi) = (Wjs72 exp 2^T~ orelse Pa ? (k) =N exp_ 2^ (10) 

in real and in the Fourier space, the form we will need later on. Taking now the form of 
the self interaction and the external fields as in Eq. ^ and El we obtain the following result 



W(r Q1 ,r 02 ) = JV^ (r i - r k ) +^^1 ( r ° 2 " rfc ) + Yl Wi + W ^ 2 (roi, r m ). (11) 

fc=l k=l k=l 
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W *= / 7^3 IPa 2 (k)| 2 n(k), (13) 



Here r-j stand for the position of the two macroins (to be distinguished from the position 
of the two centers of mass of the polymer chains roj). J-s are due to the interaction of the 
chains with external fields and can be written in the form 

^ ( r " r ') = / 70 ^ r ' r,) - (12) 

The self-interactions of the chains correspond to the terms Wj and are given by 

W) 

and finally the interactions between the two chains that can be derived in the form 

W 0l>02 (r,r') = p a? (k)n(k) Pai (-k) >. (14) 

This is all straightforward generalization of the variational theory set up previously for a 
single chain [2^. The functions Ci( r 0i) are next obtained by minimizing the upper bound 
to the exact free energy with respect to af leading to 

%(?(r 0i )N = 0£;W(r ou r O2 ). (15) 

The effective centers-of-mass free energy of the two polymer chains is finally given by the 
expression 

sinh ~ 

#Fo(roi,ro2) = 3^P log ^ - f E^i Ci Na i + PW(v 01 ,v 02 ) = 

= PT 2 (r 01 ,r 02 )+pW(r 01 ,r 02 ), (16) 

where we separated out the harmonic part of the free energy f3J- 2 (roi,ro 2 ). The first 
two terms of this variational free energy represent the entropy of the Gaussian chain 
and the last one is due to the interactions with the external fields and self-interactions. 
These are the basic equations of the Feynman - Kleinert variational theory as applied 
to the self interacting polyelectrolyte chains. They are still quite complicated because of 
the dependence on the center-of-mass coordinates roi and the final integration over these 
variables in Eq. El 

If there are no external fields that break the translational symmetry of the problem 
it can be easily seen that the dependence on ro« vanishes and the solution of the 
variational equations is straightforward. With external fields the final quite complicated 
roi integration can be obtained only numerically. In the case that ^(tqi, tq 2 ) scales with a 
positive power of iV and N is large enough, there is however an additional quite accurate 
approximation to circumvent this final integration 27 . It consists of the saddle point 
evaluation of the final integration with respect to roi, that is of an additional minimization 
of ^(roi; r 02) with respect to roi as well as ro2- This means that Eq. Elcan be written in 
an approximate form 

E (N) = [ [ d 3 r 01 d 3 r 02 e -^o(r 01 ,r 02 ) _ ^0(^02). (17) 
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where , rg 2 are given as solutions of the saddle-point condition 

dF {r* m ,r* 2 ) _ dT (r* m ,r* 02 ) _ Q 



dr oi dr 02 

Thus we obtain a simple explicit and accurate estimate for the free energy of two self- 
interacting polyelectrolyte chains in an external field, viz. T = —kT Ho(iV) ~ .Fo( r oi> r 02)- 
In what follows we will always assume that all the solutions of the variational equations 
have to be symmetric with respect to the two chains. There is no reason on the pair- 
potential level to assume otherwise. 

Since the solution of variational equations is in general quite complicated we give here 
a little preview of what exactly we will be calculating in what follows. We will show 
that grosso modo the solution of the variational problem has two branches depending on 
the relative magnitudes of the interaction with external fields and the self and mutual 
interaction of the chains. The two branches of the solution are: 

• what we call a strong coupling branch that corresponds to Cf ( r 0i) > in the varia- 
tional equation Eq. and thus to the dominance of the interactions of the chains 
with the external macroion fields, the self- and mutual interactions of the chains 
being a small perturbation. The strong coupling branch furthermore bifurcates into 
two different sub-branches depending on the solution of the additional minimization 
implied by Eq. 1181 in the weakly paired sub-branch the chains are associated each 
with its own grafting macroion and in the strongly paired sub-branch both chains 
share the two macroions on the average. 



and what we term a weak coupling branch where Ci ( r 0i) < and thus corresponds 



-2 
n 

to the case where the self and mutual interactions of the chain are dominant and the 
interactions with external macroion fields are perturbative. Coupling in both cases 
thus refers to coupling with the external macroion field. 

In both cases we can in general observe some bridging effects but they are several orders of 
magnitude stronger in the first case. Nevertheless they are always present to some extent. 
After this introductory survey of the nature of the solutions of the variational problem, 
we are ready to find these solutions explicitly. 

A note on the grafting of the chain is in order at this point. Both in the weakly paired 
as well as the strongly paired states the electrostatic adsorption energy more than the 
grafting itself determines the statistics of the chain. Grafting the chains, by fixing e.g. 
Vi(0) to be at the surface of the macroion, would change none of the conclusions reached 
below, provided of course that the size of the macroions is small compared to the size of 
the chains and that we have only one chain associated with each of the macroions. It 
would however introduce some serious complications into our formalism thus obscuring its 
straightforward interpretation. 



4 Solution of the variational equations 

We are now ready to solve the general variational equations for the case of two polyelec- 
trolyte chains with two external point macroions. As we already announced we consider 
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only symmetric solutions for which a\ = a\ = a 2 , but in general with roi 7^ i"02- This 
symmetrization will be applied to the results derived below in their final form. 

A little straihgtforward algebra then leads to the following form of the total variational 
free energy Eq. ^] 

InN 2 

pWi = -^—h(0,m) 
ira 

PW aua2 (v,r') = ^lf 1 (l\ r -T'\,Ka). (19) 

e 2 

£b = 47,-ee kT was introduced above as the Bjerrum length. Also we defined the following 
function 

f°° usinuy e~ Xu \ 
fx(y,t) = / 7~2~TW\ du = 



Jo 

IT 

4y 



e Xy2 Ue-* - e- yt Erfc(^= - tVX) - e yt Erfc(^= + tV\) \ ,(20) 



where Erfc(x) is the standard complementary error function. On the other hand the 
variational equation Eq. El can be obtained just as straihgtforwardly as, 



r 9 uri \ £ b N 



3 



2^M £ g ^|r 01 - r fc |, ^\ - 2Ng (^^, «») - Ng (0, .a) 



k=\ 

(21) 

A similar equation could be obtained also for /3^W(roi, 1*02) except that roi on the r.h.s. 
would be turned into ro2- The following new function was defined above 

g(y,t) = -^fx(y,t)\ x=1 . (22) 

What the variational equation Eq. |^ really asserts is which terms are important in 
determining the statistical conformation of the chain, i.e. af in our case. The first term 
on the r.h.s. of Eq. is due to the interactions with the macroions, the second one is 
due to the interactions between the two chains and the last one is the self interaction of 
the chains. The conformation of the chain as descibed by af is thus determined by the 
relative magnitudes of these three terms. 

The properties of the solution of the variational equation Eq. |^ first of all depend 
crucially on the sign of the r.h.s. of the Eq. 1211 thus on the fact whether £ 2 is positive or 
negative. The sign of this term tells for each value of the separation between the macroions 
whether it is the interactions with the external fields or the self-interactions of the chain 
that determine the statistical conformation of the chain. 

In view of the form of the variational ansatz Eq. |l]the positive sign corresponds to a 
general confinement of the chain if compared to the case with no interactions. We will refer 
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to the ensuing interactions between the two macroions mediated by the polyelectrolyte 
chain as strong coupling limit. In the opposite case the chain is expanded if compared to 
the case with no interactions, and we will refer to the ensuing polyelectrolyte mediated 
interactions as weaic coupling limit. Both terms will be explained further below. The final 
closure for this system of variational equations is provided by the relation between £ and 
a, Eq. El 

Since in this model the external macroions break the translational symmetry of the 
system we apply also the minimization condition Eq. ^] with respect to roi as well as 
1*02, in order to avoid the final complicated integral over the centers-of-mass of the two 
polyelectrolyte chains. This minimization introduces additional features of the solutions 
of the variational equations. Taking into account the Gaussian-like form of the function 
fl(y, t) we realize, that there are in fact two different symmetric solutions to equation Eq. 

M 

• roi = i"i and ro2 = each of the chain remains associated with its grafting 
macroion 

• roi = I'm = 2( r i+ r 2)> each chain is shared by the two macroions symmetrically. 

Here we assumed that the first chain is grafted to the first macroion while the second one 
is grafted to the second macroion. We refere to the configuration of the polyelectrolye 
chains in the first case as weakly paired and in the second case as strongly paired. The 
terms are self-explanatory: in the first two cases the chain is confined to one of the 
macroions, whereas in the second case it is confined or bound by both of them. A schematic 
representation of the solutions of the variational equations is presented on Fig. ^ 



5 Strong coupling limit 

Once again strong coupling means that external fields dominate the statistical configura- 
tion of the chain and thus £ 2 > 0. In this domain of the parameter space the effect of the 
interactions of the polyelectrolyte chain with the macroions, the term proportional to M 
in Eq. 1211 determines the overall configuration of the chain. 

The strong coupling limit entails however two different polyelectrolyte equilibrium 
states as discussed above, depending on the minimization with respect to roi, 1*02: the first 
state, stable for small values of the separation between the macroions, is due to the strong 
pairing configuration of the chain with roi = ro2 = |(ri + r2). The variational equation 
for £ in this case reads 



3 a2 



7T<2 3 



r 2 > 



v/2 



3Ng (0, ko) 



(23) 



while the corresponding free energy has the form 



(3F = 2#F 2 (|ri-r 2 | 



t B N 



ira 



2 7 / 2 M 



V2, 
2a 



ri - r 2 



67V/! (0, no) 



(24) 
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The form of the dependence ^(Iri - r 2|) is of course given implicitely via the dependence of 
£ and a. Once again the chain here is bound to both macroions and its statistical properties 
are dominated by the interaction with the charges on the macroions. One would expect 
that the polyelectrolyte mediated interactions between the macroions would be strongest 
in this case. Obviously for large enough |ri — r2| the r.h.s. of Eq. |23]can become negative, 
going first through zero. This is due to the fact that g(y,t) is a decaying function of y. At 
this point the above solution ceases to be stable and we have a transition from the strongly 
paired to weakly paired branch of the strong coupling limit. The transition depends on 
the macroion parameters such as the magnitude of their charges as well as the length of 
the chains. In this sense it represents a finite size (of the chains) effect. 

The weakly paired configuration is characterized by roi = t\ and ro2 = r 2 and is the 
stable branch at larger separations between the macroions. Here the variational equation 
for £ becomes 



3^2 _ i B N 



c 



2 ' 7ra 3 



2 3 / 2 M \g (o^+sf^ln-ral,^)] - 2Ng (^^, no) - Ng (0, no) 



(25) 

The corresponding free energy in this case can be obtained as 



(3F = #F 2 (|n-r 2 | 



£ B N 



ira 



2 5/2 M 



-Wfi \ ^ ri a T2 \ naj - 2N fx (0, m)\ . (26) 

The most important term to determine the conformation of the chain is the interaction 
with the single macroion and is thus only weakly dependent on the separation between 
them. These are the first and the last term in the r.h.s. of Eq. 1251 The separation 
dependent terms act only as a perturbation to these terms. It is thus to be expected that 
the polyelectrolyte mediated interactions will be much weaker in this case. 

Also in order to get the interaction free energy in the weakly paired configuration one 
needs in addition to subtract the terms that do not depend on the separation between 
the macroions from the total free energy. This is the standard way to get the interaction 
free energy. Again the form of the dependence J^2(| r i — r 2|) is given implicitely via the 
dependence of Q and a. 

The form of the solution of Eq. El and Eq. [25] as well as the corresponding polyelec- 
trolyte mediated interaction free energy is presented in Fig|2j We see that at small enough 
separations the chain is in the strongly paired configuration, being confined symmetrycally 
by both macroions. In this regime the external field trying to confine the chain to both 
macroions wins over the chain entropy that is expanding the chain. The entropy of the 
parts of the chain spanning the region between the macroions is quite low. Its size, as de- 
scribed by a, in this case depends monotonically on the separation between the macroions 
that are effectively stretching it. 

At the instability point, reached at a well defined value of the separation between 
the macroions, the chain entropy scores a partial victory over the interactions with the 
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macroions forcing the chain to remain close to the macroion to which it is grafted. At this 
transition the chain basically relaxes the low entropy configurations of its parts confined 
between both macroions by snapping back to macroion to which it is grafted. After that 
the size of the chain is basically determined solely through the interactions of the chain 
with its grafting macroion and remains constant with separation between the macroions. 
These conclusions reached on the basis of the two chain variational approach are very 
similar to the existing simulation data ^S] for one chain in the field of two macroions. 

This scenario of chain conformations is clearly illustrated in Figs. El and El where 
one can follow the transtition of the chain from the strongly paired to the weakly paired 
configuration via the changes of a as a function of the separation between the macroions 
for two different values of the ionic strength of the uni-valent salt solution, viz 1 mM and 
60 mM. Clearly the overall effect of the salt is to quench the magnitude of the bridging 
interaction. We note about one order of magnitude difference in the strength of the 
bridging attraction at both salt activities. The increase in salt activity also quenches the 
difference between the weakly paired and the strongly paired configuration. If we compare 
the radius oif gyration of the chain (lower graphs on Figs. |2]and|HJ) we see that at higher 
salt the separation between the macroions has a smaller effect on the size of the chain. 

The ensuing chain mediated interaction free energy follows closely the equilibrium con- 
figuration of the chain between the macroions. For a strongly paired chain the interaction 
free energy shows a pronounced attractive contribution stemming from the coupling MxN 
term in Eq. 1241 At the point where the chain snaps from the strongly paired to the weakly 
paired configuration there is also a corresponding jump in the free energy due to much 
less pronounced chain mediated interactions. It is interesting to analyze the asymptotic 
behavior of the interaction free energy Eq. in the weakly paired branch of the solu- 
tion. Expanding t) for large values of y we see that all the chain dependent parts 
of the free energy finally just renormalize the direct Debye-Hiickel interactions between 
the macroions . Thus asymptotically instead of a Debye-Hiickel interaction of strength 
proportional to M 2 we simply end up with its strength being proportional to (M — N) 2 . 
The chain snapped back to the grafting macroion and simply renormalized its charge. 

There is one interesting remark that we can make here. Gurovitch and Sans 28 studied 
polyelectrolyte adsorption of a single chain to a (point) charged macoion. Their case thus 
correspond to an weakly paired branch at infinite separation between the macroions in 
our terminology, that would correspond to Eq, [23 with |ri — rg| — ► oo leading to 



polyelectrolyte chain can adsorb until its charge (or the number of monomers) becomes 
equal to N = 2 3//2 M, there is thus maximal overcharging in the amount of 2 3 / 2 « 2.83 
which is indeed very close to the value derived by them by a completely different method 
viz. 15/6 = 2.5. Though the approach of these authors has been criticized [2Hj our results 
are more then consistent with theirs. More could be said on polyelectrolyte adsorption and 
overcharging [2S] but we will focus here strictly on the interaction i.e. bridging aspects of 
the problem. 




(27) 



Clearly in the vanishing salt limit kcl 



0, which is in fact the case treated in [2B]j the 
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6 Weak coupling limit 



In this case the effect of electrostatic self-interaction of the chain, the term proportional to 
N in Eq. I18| determines the overall configuration of the chain. If the effect of the external 
fields would be indeed negligible we have shown in a previous publication |21| that the 
electrostatic interactions would stiffen up the chain and give it a rod-like appearance 
quantified by the scaling a ~ N. We expect that even with external fields originating at 
the macroions the chain will essentially assume this type of extended configurations in this 
limit, modified by the perturbative effect of both macroions. Simulations of single chain 
adsorption |23| are completely consistent with this picture since for large N protruding 
rodlike tails are observed that correspond to electrostatically stiffened portions of the 
chain. Calculations of Nguyen and Shklovskii j2H] also lead to the same qualitative picture 
of chain adsorption in this limit. Again here we are not interested in adsorption per se 
but rather in the bridging effects in interaction between the macroions, so we will skip the 
detailed comparison of our work with polyelectrolyte adsorption studies. 

With the external fields present the stiffening of the chain depends on the details i.e. 
strength of the term proportional to M in the variational equation Eq. El an d is difficult to 
quantify in scaling terms, unless the effect of the macroions is thoroughly negligible, which 
is not the situation we are trying to investigate. Since here the effect of the macroions is 
small the polyelectrolyte chains can never be strongly paired by both macroions. We thus 
remain solely with the weakly paired configuration of the chains due to the grafting to the 
macroions. The solution of the variational equation Eq. ^] thus only has one branch in 
this case and is given by 



3 (-2 _ e B N 

2^ 



Wg (^d, «a) + Ng (0, «a) - 2 3 / 2 M [g (o, + g (^|n - r 2 |, ^ 

'(28) 

Clearly this equation is obtained by the substitution £ — > iQ from Eq. 1181 This transfor- 
mation should be taken into account also in Eq. |H1 leading to 



.2 



1 (Ct 2 N 



a 



where now £{x) = - — cotx. The interaction with the macroion, the M term in the Eq. 
1281 can modify the value of the size of the chain but it has no effect any more on the 
stability of the solution. The numerical solutions to Eq. |25]are presented in Figs. 2]and[SJ 
Clearly the size of the polyelectrolyte chain in this case shows no discontinuities though it 
is still, to a lesser extent than before, effected by the positions of the two macroions. The 
term weak coupling thus seems appropriate for the behavior of the chain in this region of 
the parameter space. 

The free energy is now given by an equation similar to Eq. [2^1 but with the change 
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i( well taken into account in ^(l 1 "! — r 2|)- It leads to the following result 



6 log 




+ 3 C 2 iVa 2 




(30) 



Again because the solution of the variational equations here remains on a single branch 
all the time, showing no jump from one stable branch to another one ("snapping" of the 
chain). The free energy shows no discontinuities either, though it still depends on the 
separation between macroions. 

This state of affairs introduces new features in the polyelectrolyte mediated interac- 
tions. First of all there is a clearly discernible, see Fig. |1J repulsion at smaller separations. 
It is due to the interpenetrating "coronas", i.e. extended configurations of the grafted 
chains, on approach of the macroions. If there would be many chains grafted to both 
macroions this incipient repulsion would develop into a full-blown brush repulsion regime. 
Since we have only one chain per macroion the effect of "coronal" interpenetration is rather 
weak, but nevertheless clearly discernible. Its range depends on the size of the chain N 
as well as the amount of salt which regulates the overall extension of the chain. Fig. 
clearly shows that salt quenches the "coronal" repulsion. 

It is only at larger separations, see the inset of Fig. |1J that residual, indeed very 
weak coupling, is finally discerned. The electrostatically extended chains can still make 
weak bridges to the other macroions but since this can happen only at sufficiently large 
separations, the ensuing bridging is much attenuated. On adding the salt this effect is 
displaced towards smaller separations because the extent of the chain is deminished by 
the salt as well. For larger ionic strengths we are thus left with weak coupling at smaller 
separations, see Fig. where it is nevertheless stronger than in small salt, compare again 
inset to Fig. ^ Bridging, no matter what its source, is of course stongly dependent on 
the separations between the macroions and is in general stronger for smaller separations. 
Clearly the weak coupling interaction in the case of large salt resembles much more " sticky 
macroions" then (relatively) long range bridging interaction. 

In the weak coupling limit there is thus an additional feature stemming from the 
polyelectrolyte mediated interactions which is due to "coronal" interpenetration and marks 
the incipient brush repulsion that would be developed fully if more chains were grafted 
to each macroions. The (weak) bridging attraction in this case is overall small and is in 
constant competition with "coronal" interpenetration interactions. 

7 Discussion 

The polyelectrolyte bridging interaction analyzed here is obviously very rich in its features 
and depends crucially on the region of the parameter space under investigation. We showed 
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that attractive bridging interactions effectively come in two varieties: the strong coupling, 
where interactions between the chains and the macroions are dominant, and the weak 
coupling, where self-interactions of the chains are dominant. Bridging interaction can be 
obtained in both cases but is a couple of orders of magnitude larger in the strong coupling 
limit. It the weak coupling limit the attractive bridging interactionis are overwhelmed by 
the incipient electrostatic repulsions between the overlapping polyelectrolyte "coronas" as 
well. In both cases however the effective polyelectrolyte mediated interaction potential is 
strongly non-monotonic and shows pronounced variation with respect to the length of the 
polyelectrolyte chains and the screening length of the underlying electrostatic interactions. 

The variational formulation of the bridging interaction problem which lies at the basis 
of our approach, has several advantages as well as drawbacks. The main feature of the 
formalism is that it allows for the transition between the strong pairing and weak pairing 
states of the polyelectrolyte chains which is clearly a finite size effect and would thus be 
missed on the ground-state dominance level. The latter has been used succesfully for the 
polyelectrolyte mediated interactions between macroscopic surfaces with intervening long 
polyelectrolyte chains |13| I16j . Another feature of the variational formulation is that it 
can describe the snapping of the chain at large enough macroion separations which is the 
most important feature of a finite chain size effect. We need to reiterate again that this 
phenomenon is absent in the ground state dominance ansatz. The snapping of the chain 
from the configuration where it is partitioned between both macroions to one, where it 
is adsorbed to a single one, is of course due to an interplay and balance between chain 
adsorption energy, chain self-interaction energy and configurational entropy of the chain. 
The balance depends on the size of the chain and the separation between the macroions 
and leads to an abrupt transition between the two configurations that has a well discerned 
imprint also on the polyelectrolyte mediated interaction between the macroions. In general 
one sees bridging only for chain cofigurations where it is partitioned by the two macroions 
symmetricaly, i.e. in what we dubbed the strong pairing configuration. 

The attractive bridging interaction is typicaly about a hundred times stronger in the 
strong coupling limit if compared to its weak coupling counterpart (whence the designation 
of the two limiting cases). This is intuitively quite easy to grasp, since one can expect 
strong coupling to emerge only when the interactions between the chain and the macroions 
dominate the statistical properties of the system. If however the dominant interactions 
in the system are self-interactions of the chains, the polyelectrolytes clearly mediate only 
insignificantly the interactions between the macroions. In this case a much more important 
feature of the interaction is the interpenetration of the polyelectrolyte "coronas", that can 
in some cases lead to pronounced repulsions between "dressed" macroions. These repusions 
would clearly stabilize the macroion interactions. 

The main drawback of the present analysis of the bridging interaction problem, apart 
from it being a purely pairwise additive formulation, is the linearized (Debye-Hiickel) form 
of electrostatics. All non- linear effects are thus a priori excluded. On this level the main 
effect of the salt is to attenuate the bridging interaction as well as the repulsive interaction 
between the polyelectrolyte "coronas". In this respect the variational approach is inferior 
to the ground-state dominance ansatz. One possible way out would be to formulate also 
the electrostatic part of the problem on a variational level where all the Debye-Hiickel 
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parameters would be determined self-consistently. We leave this exercise for future work. 

Another importan omission of our method is the size of the macroions that does not 
feature explicitly in our formulation. The finite size of the chain in weakly paired or 
strongly paired configuration clearly showed by our numerical results (see the lower graphs 
in Figs. |2]and|3J) is thus not due to the finite size of the adsorbing macroions as in more 
realistic simulations \2'A\ but is an entropy-energy competiton effect: high adsorption 
energy vs. low configurational entropy in the weakly paired state, leading to the finite 
size of the weakly paired state even with a point adsorbing macroion. The omission of the 
finite size of the macroion tends to overestimate the polyelectrolyte mediated interactions 
and underestimate the size of the chain in the weakly paired as well as strongly paired 
configurations. This type of finite size effects can be straightforwardly incorporated into 
the statistics of free, non-interacting chains |30| . but would be unfortunately difficult to 
incorporate into the Feynman-Kleinert variational method and were thus ignored in our 
formulation. Alternative approaches would thus have to be considered |31j . 

At this stage it does not seem reasonable to compare the second virial coefficient de- 
rived from our calculation of the effective pair interactions with the experiment on NCPs 
[TU] . There would be just too many adjustments that one would have to put in by hand 
but that would have a crucial effect on the ensuing numerical results: the effective charge 
of the macroion due to non-linear Poisson-Boltzmann effect, the effective charge of the 
chains which depends strongly on the local ionic equilibrium of the dissociable amino acid 
groups, the effective length of the chains that are free enough to behave as flexible poly- 
electrolytes etc. Nevertheless, if one chooses to ignore all these additional complications 
the virial coefficient in the strong coupling limit comes out always monotonically (decreas- 
ing) dependent on the screening length, i.e. ionic strength. No non-monotonic effects, 
of the type that feature so prominently in experimetal results |10| . are ever seen for any 
reasonable values of parameters (i.e. the charges on the macroions, the length of the chain, 
the charges on the chain). Basing our conslusion on the anlysis presented above, we are 
inclined to believe that the polyelectrolyte bridging itself never leads to non-monotonic 
second virial coefficient. One nevertheless has to keep in mind, that our model calculation 
is based on many constraints that are not entirely realistic. 

Our work represents an alternative formulation of the polyelectrolyte bridging inter- 
action between two small macroions, on the two particle level. The finite size effects 
of the chain length make obviously a strong imprint on the bridging interaction. These 
effects have not been studied previously analytically and are missed by the more pop- 
ular ground-state dominance mena-field approach. In this respect we believe our work 
adds an important feature to our understanding of the phenomenon of polymer mediated 
interactions. 
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Figure captions 
Figure [I] 

A schematic representation of the model and the nature of the variational solutions. 
L.h.s.: the strong coupling solution branch. The conformation of the chain is determined 
mostly by the interaction with the two macroions. It can be either in the weakly (upper) 
or strongly paired (lower drawing) configuration, depending on the separations between 
the macroions. The attractive bridging interaction is strongest in the strongly paired 
configuration. R.h.s.: the weaic coupling solution branch. The conformation of the chain 
is determined mostly by the self-interactions of the chains. They are in the weakly paired 
configuration for any separation between the macroions. On approach the interpenetration 
of the monomer clouds (the "coronas") leads to prevalently repulsive interactions. Note 
that the macroions are treated as point aprticles in the claculation and that their size is 
thus exagerated in this drawing. 

Figure El 

(3J-q and a in the strong coupling limit. The upper graph shows the dependence of (3J-Q 
for M = 100 for N = 10 (losenge), 30 (square), 100 (circle) at ImM on the separation 
between the macroions |ri — r2|. The lower graph shows the dependence of the size of the 
chain a on the separation between the macroions for the same values of parameters. The 
length of the chain N obviously determines the separation between macroions where the 
snapping of the chain between the strongly paired and weakly paired states occurs. In 
the strongly paired configuration we have a well developed regime of attractive bridging 
interactions leading to a an effective screened coulomb repulsion in the weakly paired 
regime. In all cases £ = 10. The bold line represents the pure Debye-Hiickel interactions 
between the macroions. 

Figure |3] 

(3J-q and a in the strong coupling limit. The upper graph shows the dependence of 
/3J-q for M = 100 for N = 10, 30, 100 at 60mM on the separation between the macroions 
| ri — r2 1 . The lower graph shows the dependence of the size of the chain a on the separation 
between the macroions for the same values of parameters. The bold line represents the 
pure Debye-Hiickel interactions between the macroions. 

Figure |3] 

PJ-q and a in the weak coupling limit for M = 40 and iV = 50, 20, 10 at ImM. There 
exists only an weakly paired state in this case and the effect of the external fields of the 
macroions is much less pronounced then in the strong coupling limit. The length of the 
chain N determines primarily the steric repulsion effect due to the interpenetration of stiff- 
ened grafted chains, the "coronas" of both macroions, on close approach. The inset shows 
the residual very weak coupling interaction at large separations between macroions with 
N = 100. The size of the chain tends to grow slightly on approach of the macroions be- 
cause of the interpenetration of the "coronas" of both macroions. The bold line represents 
the pure Debye-Hiickel interactions between the macroions. 

Figure |S] 
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(5Tq and a in the weak coupling limit for M = 40 and N = 50, 20, 10 at 60mM. Clearly 
again the overall effect of the salt is to quench the effects of the macroion fields. There 
is weak residual bridging at very small separation now, note the scale of the interactions 
energy, that goes smoothly into a weak "corona" repulsion at larger separations. The 
details of both the residual bridging as well as "corona" repulsions depend on the values 
of M and N. The bold line represents the pure Debye-Hiickel interactions between the 
macroions. 
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Figure 1: A schematic representation of the model and the nature of the variational so- 
lutions. L.h.s.: the strong coupling solution branch. The conformation of the chain is 
determined mostly by the interaction with the two macroions. It can be either in the 
weakly (upper) or strongly paired (lower drawing) configuration, depending on the sep- 
arations between the macroions. The attractive bridging interaction is strongest in the 
strongly paired configuration. R.h.s.: the weaic coupling solution branch. The conforma- 
tion of the chain is determined mostly by the self-interactions of the chains. They are in 
the weakly paired configuration for any separation between the macroions. On approach 
the interpenetration of the monomer clouds (the " coronas" ) leads to prevalently repulsive 
interactions. Note that the macroions are treated as point aprticles in the claculation and 
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Figure 2: (3J-q and a in the strong coupling limit. The upper graph shows the dependence 
of PFo for M = 100 for N = 10 (losenge), 30 (square), 100 (circle) at ImM on the 
separation between the macroions |ri — r2|. The lower graph shows the dependence of 
the size of the chain a on the separation between the macroions for the same values 
of parameters. The length of the chain N obviously determines the separation between 
macroions where the snapping of the chain between the strongly paired and weakly paired 
states occurs. In the strongly paired configuration we have a well developed regime of 
attractive bridging interactions leading to a an effective screened coulomb repulsion in the 
weakly paired regime. In all cases £ = 10. The bold line represents the pure Debye-Hiickel 
interactions between the macroions. 
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Figure 3: (3J-q and a in the strong coupling limit. The upper graph shows the dependence 
of (3J-q for M = 100 for N = 10, 30, 100 at 60mM on the separation between the macroions 
| r i — T2 1 . The lower graph shows the dependence of the size of the chain a on the separation 
between the macroions for the same values of parameters. The bold line represents the 
pure Debye-Hiickel interactions between the macroions. 
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Figure 4: /3jF and a in the weak coupling limit for M = 40 and iV = 50, 20, 10 at ImM. 
There exists only an weakly paired state in this case and the effect of the external fields 
of the macroions is much less pronounced then in the strong coupling limit. The length 
of the chain N determines primarily the steric repulsion effect due to the interpenetration 
of stiffened grafted chains, the "coronas" of both macroions, on close approach. The inset 
shows the residual very weak coupling interaction at large separations between macroions 
with N = 100. The size of the chain tends to grow slightly on approach of the macroions 
because of the interpenetration of the "coronas" of both macroions. The bold line repre- 
sents the pure Debye-Hiickel interactions between the macroions. 
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Figure 5: f3To and a in the weak coupling limit for M = 40 and TV = 50, 20, 10 at 60mM. 
Clearly again the overall effect of the salt is to quench the effects of the macroion fields. 
There is weak residual bridging at very small separation now, note the scale of the inter- 
actions energy, that goes smoothly into a weak "corona" repulsion at larger separations. 
The details of both the residual bridging as well as "corona" repulsions depend on the 
values of M and N. The bold line represents the pure Debye-Huckel interactions between 
the macroions. 
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